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Abstract 

In this paper we study minimum cut and maximum flow problems on planar graphs, both 
in static and in dynamic settings. First, we present an algorithm that given an undirected 
planar graph computes the minimum cut between any two given vertices in 0(nloglogn,) time. 
Second, we show how to achieve the same O(nloglogn) bound for the problem of computing 
maximum flows in undirected planar graphs. To the best of our knowledge, these are the first 
algorithms for those two problems that break the 0(n log n) barrier, which has been standing 
for more than 25 years. Third, we present a fully dynamic algorithm that is able to maintain 
information about minimum cuts and maximum flows in a plane graph (i.e., a planar graph 
with a fixed embedding): our algorithm is able to insert edges, delete edges and answer min-cut 
and max- flow queries between any pair of vertices in 0(n 2//3 log 3 n) time per operation. This 
result is based on a new dynamic shortest path algorithm for planar graphs which may be of 
independent interest. We remark that this is the first known non-trivial algorithm for min-cut 
and max-flow problems in a dynamic setting. 
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1 Introduction 



Minimum cut and maximum flow problems have been at the heart of algorithmic research on graphs 
for over 50 years. Particular attention has been given to solving those problems on planar graphs, 
not only because they often admit faster algorithms than general graphs but also since planar 
graphs arise naturally in many applications. The pioneering work of Ford and Fulkerson [3l 0], 
which introduced the max-flow min-cut theorem, also contained an elegant algorithm for comput- 
ing maximum flows in (s, t)-planar graphs (i.e., planar graphs where both the source s and the sink 
t lie on the same face). The algorithm was implemented to work in 0{n log n) time by Itai and 
Shiloach [9]. Later, a simpler algorithm for the same problem was given by Hassin [6], who reduced 
the problem to single-source shortest path computations in the dual graph. The time required 
to compute single-source shortest paths in (s,t)-planar graphs was shown to be 0{n\J\ogn) by 
Frederickson [5] and later improved to 0(n) by Henzinger et al. [8]. As a result, minimum cuts and 
maximum flows can be found in 0(n) time in (s,t)-planar graphs. 

Itai and Shiloach [9] generalized their approach to the case of general planar (i.e., not only 
(s, t)-planar) graphs, by observing that the minimum cut separating vertices s and t in a planar 
graph G is related to the minimum cost cycle that separates faces f s and ft (corresponding to 
vertices s and t) in the dual graph. The resulting algorithm makes 0{n) calls to their original 
algorithm for (s,t)-planar graphs and thus runs in a total of 0(n 2 log n) time. In the case of undi- 
rected planar graphs, Reif [13] improved this bound by describing how to find the minimum cost 
separating cycle with a divide-and-conquer approach using only O(logn) runs of the (s,t)-planar 
algorithm: this yields an 0(n log 2 n) time to compute a minimum cut for undirected planar graphs. 
Later on, Frederickson [5] improved the running time of Reif's algorithm to O(ralogn). The same 
result can be obtained by using more recent planar shortest path algorithms (see e.g., [8]). Hassin 
and Johnson [7J extended the minimum cut algorithm of Reif to compute a maximum flow in only 
0{n log n) additional time: this implies an undirected planar maximum flow algorithm that runs 
in O(ralogre) time as well. In summary, the best bound known for computing minimum cuts and 
maximum flows in planar undirected graphs is O(ralogn). 

The first contribution of this paper is to improve to 0(n log log n) the time for computing 
minimum cuts in planar undirected graphs. To achieve this bound, we improve Reif's classical 
approach [13] with several novel ideas. To compute a minimum s-t cut in a planar graph G, we first 
identify a path tt between face f s (corresponding to vertex s) and face ft (corresponding to vertex 
t) in the dual graph Gd- Next, we compute a new graph G n as follows. We cut Gd along path tt, 
so that the vertices and edges of tt are duplicated and lie on the boundary of a new face: another 
copy of the same cut graph is embedded inside this face. We show that minimum separating cycles 
in correspond to some kind of shortest paths in G n . Applying a divide-and-conquer approach 
on the path tt yields the same 0{n log n) time bound as previously known algorithms [7J 113] . 
However, our novel approach has the main advantage that it allows the use of any path tt in the 
dual graph Go, while previous algorithms were constrained to choose tt as a shortest path. We 
will exploit the freedom implicit in the choice of this path tt to produce a faster 0(n log log n) time 
algorithm, by using a suitably defined cluster decomposition of a planar graph, combined with 
the Dijkstra-like shortest path algorithm by Fakcharoenphol and Rao [2]. Our second contribu- 
tion is to show that also maximum flows can be computed in undirected planar graphs within the 
same 0(n log log n) time bound. We remark that this is not an immediate consequence of our new 
minimum cut algorithm: indeed the approach of Hassin and Johnson [7j to extend minimum cut 
algorithms to the problem of computing maximum flows has a higher overhead of O(ralogra). To 
get improved maximum flow algorithms, we have to appropriately modify the original technique of 
Hassin and Johnson [7J. To the best of our knowledge, the algorithms presented in this paper are 
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the first algorithms that break the 0(n log n) long-standing barrier for minimum cut and maximum 
flow problems in undirected planar graphs. 

As our third contibution, we present a fully dynamic algorithm that is able to maintain infor- 
mation about minimum cuts and maximum flows in a plane graph (i.e., a planar graph with a fixed 
embedding): our algorithm is able to insert edges, delete edges and answer min-cut and max- flow 
queries between any pair of vertices in 0(n 2 / 3 log 3 n) time per operation. This result is based on 
the techniques developed in this paper for the static minimum cut algorithm and on a new dynamic 
shortest path algorithm for planar graphs which may be of independent interest. We remark that 
this is the first known non-trivial algorithm for min-cut and max-flow problems in a dynamic setting. 

2 Minimum Cuts in Planar Graphs 

Let G = (V, E, c) be a planar undirected graph where V is the vertex set, E is the edge set and 
c : E — >• 1Z + is the edge capacity function. Let the planar graph G be given with a certain embed- 
ding. Using the topological incidence relationship between edges and faces of G, one can define the 
dual graph Gd = {F, Ed, cb) as follows. Each face of G gives rise to a vertex in F. Dual vertices /i 
and f2 are connected by a dual undirected edge eo whenever primal edge e is adjacent to the faces 
of G corresponding to f\ and f2- The weight cd(&d) of the dual edge &b is equal to the weight 
c(e) of the primal edge: co(eo) is referred to as the length of edge en- In other terms, the length 
of the dual edge is equal to the capacity of the primal edge e. In the following, we refer to G as 
the primal graph and to Gd as its dual. Throughout the paper we will refer to vertices of the dual 
graph Gd interchangeably as (dual) vertices or faces. Note that Go can be embedded in the plane 
by placing each dual vertex inside the corresponding face of G, and placing dual edges so that each 
one crosses only its corresponding primal edge. Thus, the dual graph is planar as well, although 
it might contain multiple edges and self loops. In simpler terms, the dual graph Go of a planar 
embedded graph G is obtained by exchanging the roles of faces and vertices, and G and Gd are 
each other's dual. Figure [1] shows an embedded planar graph G and its dual Gd- 

Let s and t be any two vertices of G (not necessarily on the same face). We consider the 
problem of finding a minimum cut in G between vertices s and t. Let C be a cycle of graph G: we 
define the interior of C, denoted by int(C), to be the region inside C and including C in the planar 
embedding of the graph G. We can define the exterior ext(C) of the cycle C in a similar fashion. A 
cycle C in Gd is said to be a cut-cycle if int(C) contains exactly s but not t. The following lemma 
was proven by Johnson [10j . 

Lemma 1. A minimum s-t cut in G has the same cost as a minimum cost cut-cycle of Gd- 

The lemma follows by the observation that for any cut-cycle C the faces of Gd inside int(C) 
give a set of vertices S in G which defines a cut separating s and t. Note that Lemma [TJ gives an 
equivalence between min-cuts in the primal graph G and minimum cost cut-cycles in the dual graph 
Gd- By using a divide-and-conquer approach, this equivalence can be turned into an efficient al- 
gorithm for finding flows in undirected planar graphs |13| . The resulting algorithm, combined with 
more recent results on shortest paths in planar graphs [8], is able to work in a total of 0(n log n) 
time. However, this approach seems to inherently require 0(n log n) time, and does not seem to 
leave margin for improvements. In the next section, we will present a completely different and more 
flexible approach, which will yield faster running times. 

2.1 Computing Min-Cuts 

Let f s and ft be arbitrary inner faces incident to s and to t respectively. Find any simple path ir 
from f s to ft in Gd- The path ir can be viewed as connecting special vertices in the dual graph 
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corresponding to s and t. Hence, any s-t cut needs to cross this path, because it splits s from t. 
Let ir traverse dual vertices f\, . . . , fk, where f\ = f s and fk = ft- Let us look at the path ir as a 
horizontal line, with f s on the left and ft on the right (see Figure (2{b)). An edge ep ir in Gd such 
that &d is incident to some face fi, 1 < i < k, can be viewed as connected to fi from below or from 
above. We now define a new graph G' n , by cutting Gd along path ir, so that the vertices and edges of 
7T are duplicated and lie on the boundary of a new face. This is done as follows. Let ir' be a copy of 
ir, traversing new vertices /[, f 2 , . . . , f' k . Then is the graph obtained from Gd by reconnecting to 
edges entering fi from above, 1 < i < k (see Figure Efc)). Let G" be a copy of G' n . Turn over the 
graph G" to make the face defined by ir and ir' to be the outer face (see Figure [2(d)). Now, identify 
vertices on the path ir' (respectively ir) in G" with the vertices on the path ir (respectively ir') in 
G' n . Denote the resulting graph as (see Figure Efe)). Note that the obtained graph G w is planar. 

We define an fa-cut-cycle to be a cut-cycle in Gd that includes face fi and does not include any 
face fj for j > i. The proof of the following lemma is immediate. 

Lemma 2. Let Ci be a minimum fi-cut- cycle in Gd for i = 1, . . . , k. Then Ci with minimum cost 
is a minimum cut-cycle in Gd- 

A path p between fi and f[ in G n is said to be fi- separating if p contains neither fj nor /j, for 
j > i. We say that a cycle C in Gd touches path ir in face fi if two edges on C incident to fi go 
both up or both down, whereas we say that C crosses ir in face fi if one of these edges goes up, 
whereas the other goes down. 

Lemma 3. The cost of a minimum fi-cut- cycle in Gd is equal to the length of a shortest fi- 
separating path in G v . 

Proof. Let C be some /j-cut-cycle in Gd'- we show that there must be some /j-separating path p in 
G-n having the same cost as C. Note that the /j-cut-cycle C must either cross or touch the path ir 
in face fi. First, assume that C crosses irm.fi. Note that in this case C has to cross the path ir an 
even number of times (excluding /j), as otherwise C would not separate s from t (see Figure El^a)). 
We can go along C starting from fi in the graph G n and each time when C crosses ir in Gd, we 
switch between G' n and G" in G^. Hence, due to parity of the number of crossings with ir, this will 
produce a resulting path p in G w which must end in f[ (see Figure E{b)). Second, assume that C 
touches 7r in fi. Then C has to cross path ir an odd number of times (see Figure E|c)). Again if we 
trace C starting from /, in G n we will produce a path p in G n ending up in f[ (see Figure E(d)). 
Moreover, since C is a /j-cut-cycle in Gd, it cannot contain by definition any face fj, for j > i: 
consequently, in either case the resulting path p in G n will contain neither fj nor f'-, for j > i. 

Conversely, let p be some /j-separating path in G^: we show that there must be some /j-cut- 
cycle C in Gd having the same cost as p. First, assume that p enters /j and f[ using edges from the 
same graph, i.e., either G' v or G". If this is the case, then the two edges must be such that one is 
from above and the other is from below. By tracing p in Gd we obtain a cycle C that crosses ir in 
fi. This cycle does not need to be simple, but it has to cross ir an odd number of times (including 
the crossing at fi), as each time when C crosses ir the path p must switch between G'^ and G". 
Hence, C must be a cut-cycle. Second, assume that p enters /j and f[ using edges from different 
graphs, i.e., one edge from G' n and the other edge from G". If this is the case, then the two edges 
must be both from above or both from below. Consider first the case where p leaves fi from above 
using an edge from G", and enters f[ from above using an edge from G' n (see Figure E^d)). This 
time we get a cut-cycle C touching ir in /j and intersecting ir an odd number of times. The case 
where p leaves /j from below using an edge from G' v , and enters f[ from below using an edge from 
G" is completely analogous (see Figure 0(f)). Once again, since p is an /j-separating path in G n , 
in any case the resulting cut cycle C will not contain fj, for j > i. □ 
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Note that, having constructed the graph G n , we can find the shortest /j-separating path by 
simply running Dijkstra's algorithm on G n where we remove all faces fj and /j for j > i. By Lem- 
mas [Tj [2] and [3l the linear-time implementation of Dijkstra's algorithm known for planar graphs [8] 
implies an 0(n 2 ) time algorithm for computing minimum cuts in planar graphs. There is a more 
efficient way of computing minimum cuts that goes along the lines of Reif's recursive algorithm j!3j . 
Before describing this approach, we need to prove some non-crossing properties of /j-cut-cycles. 

Lemma 4. For i < j, let Ci be a minimum /j- cut- cycle in Go, end let Cj be a minimum f j-cut- 
cycle in Go- Then there exists a cut cycle C C int{Ci) n int(Cj) in Go such that c(C) < c(Cj). 

Proof. If Ci C int(Cj) then C = Ci and the lemma follows trivially. On the other hand, it is 
impossible that Cj C int(Ci) because in this case Cj would include some face /& on tt for k > j > i. 

The only possibility left is that Ci <£. int(Cj) and Cj ^ int(Ci). In this case, there must exist 
a subpath of Ci from f a to ft, such that pi intersects int{Cj) only at f a and Let pj be the 
subpath of Cj going from /„ to fb (see Figure H]). We claim that c{pj) < c{pi). Indeed, suppose by 
contradiction that c(pj) > c(pi), and let Cj be the cycle obtained from Cj after replacing path pj 
with pi. Then the cycle Cj is shorter than Cj. Moreover, since Cj does not include any face fk for k > 
i, also Cj cannot include any face fk for k > j > i. As a consequence, the cycle Cj is an /j-cut-cycle 
in Go, with c(Cj) < c(Cj), contradicting our assumption that Cj is a minimum /j-cut-cycle in Gd- 

Since c{pj ) < c(pi), replacing path pi on Cj with the path pj yields a cycle C[, with c{C[) < c(Cj) . 
As long as C| ^ int(Cj) and Cj ^ int(Cl) we can repeat this procedure. At the end, we will obtain 
a cycle C C int(Ci) D int(Cj), such that c(C) < c(Cj). Moreover, the obtained cycle C will be a 
cut-cycle, as int(Ci) n int{Cj) contains s. □ 

Lemma 5. For i < j, Zei Cj be a minimum fi-cut-cycle in Go, and let Cj be a minimum f j-cut- 
cycle in Go- Then, for some i' < i, there exists a minimum -cut-cycle Cj' C int(Cj) such that 
c(Ci,)<c(Ci). 

Proof. Consider the cut-cycle C contained inside int{Ci) n int{Cj) as given by LemmaHl for which 
we know that c(C) < c(Cj). As C C int(Ci), C cannot contain any fk such that k > i. Hence, it is 
an /j/-cut-cycle for some i 1 < i. The minimum /j'-cut-cycle Cj' is shorter than C and consequently 
it is shorter than Cj. If Cj/ C int(Cj), the lemma follows. Otherwise, we can apply Lemma 2] to 
produce another minimum /j//-cut-cycle contained in int(Cj) such that c(Cj") < c(Cj), for some 
i" <»'<». □ 

The above lemma shows that each computed cut splits the graph into two parts, the interior 
and the exterior part of the cut, which can be handled separately. Hence, we can use a divide 
and conquer approach on the path n. We first find a minimum cut-cycle that contains the middle 
vertex on the path it. Then we recurse on both parts of the path, so there will be O(logra) levels of 
recursion in total. In this recursion we will compute minimum cost s-t cuts for all vertices on tt. By 
Lemma O we compute a minimum /j-cut-cycle for some /j by finding a shortest path in the planar 
graph G^. Then we need to divide the graph into the inside and the outside of the cut. The vertices 
on the minimum cut need to be included into both parts, so we need to take care that the total size 
of the parts does not increase too much. This can be done in a standard way as described by Reif |13| 
or by Hassin and Johnson jTj. Their technique guarantees that on each level of the recursion the 
total size of the parts is bounded by 0(n). Hence, using the 0(n)-time algorithm [8] for shortest 
paths we get an 0(n log n) time algorithm for finding minimum cuts in planar graphs. Although 
this approach yields the same time bounds as previously known algorithms [3[T3], it has the main 
advantage to allow the use of any path tt in Go, while previous algorithms were constrained to 
choose tt as a shortest path. As shown in the next section, we can even allow the path tt to be 



4 



non-simple. In Section [3] we will show how to exploit the freedom implicit in the choice of the path 
tt to produce a faster 0(n log log n) time algorithm for finding minimum cuts in planar graphs. 

2.2 Using Non-simple Paths 

Let 7T = (vi, . . . , v#) be a non-simple path and let v be a vertex appearing at least twice on the path 
7r, i.e., v = v i = Vj for some i < j. We say that the path tt is self-crossing in v if the edges incident 
to v on tt appear in the following circular order (vi-i,Vi), (vj-\,Vj), (1^,1^+1), (vj,Vj+i) in the 
embedding. We say that a path tt is self-crossing if tt is self-crossing in some vertex v. Otherwise 
we say that tt is non-crossing. If a vertex v appears at least twice on a non-crossing path tt than we 
say that tt touches itself in v. In the previous section we assumed that the path tt from / s to ft in 
Go was simple, now we will show that we only need the weaker assumption that tt is non-crossing. 

In order to work with non-crossing paths we will modify the graph Gpj to make the path tt 
simple. Let v = Vi be a face where tt touches itself in Go- Note that there is no other edge from tt 
between the edges (vf-i, Vi), (vi,Vi + i) in the circular order around v given by the embedding. Take 
all edges E v incident to v that are embedded between and including the edges (vi-±,Vi), «t+i)- 
Now, add a new face v' to Gjj and make the edges E v to be incident with v' instead of v. Finally 
connect v with v' using an undirected edge of length zero (see Figure [5]). Let tt be a non-crossing 
path: perform this vertex-splitting operation until tt becomes a simple path. This produced a new 
graph Go,n- Note that this transformation does not change the lengths of cut-cycles, so we get the 
following observation. 

Corollary 6. The lengths of minimum fi-cut-cycles in Go and Go,n Q- re the same. 

As a result, if tt is a non-simple non-crossing path, instead of running our algorithm on Gjj, we 
can compute the graph Gd,-k and run our algorithm on Gd,-k- 

3 Cluster Partitions 

Our algorithms are based on a particular cluster decomposition. We start by presenting the ideas 
behind this cluster decomposition and then show how this decomposition can be effectively ex- 
ploited for computing min-cuts. From now on, we assume that we are given a graph for which 
the dual graph has vertex degree at most three. This is without loss of generality, since it can be 
obtained by triangulating the primal graph with zero capacity edges. 

Let n be the number of vertices of G. We first define a cluster partition of G into edge clusters 
which will be used by our algorithm. In the cluster partition the graph is decomposed into a set 
of cluster V such that each edge of G is included into exactly one cluster. Each cluster contains 
two types of vertices: internal vertices, and border vertices. An internal vertex is adjacent only to 
vertices in the same cluster, while a border vertex is adjacent to vertices in different clusters. A 
hole in a cluster is any face (including the external face) containing only boundary vertices. We 
denote by dP the set of border vertices of a cluster P. We define an r ■-partition of an n-node graph 
to be a cluster partition that contains at most 0(j) clusters, each containing at most r vertices, 
0(y/r) border vertices and a constant number of holes. The proof of the following lemma is in 
Appendix |Al It is based on ideas of Frederickson [5], who constructed a similar partition without 
the bound on the number of holes. 

Lemma 7. An r -partition of an n-node planar graph can be computed in 0(n log r + log n) time. 

We use the r-partition to define a representation for shortest paths in a graph that has similar 
number of edges, but fewer vertices. In order to achieve this, we use the notion of dense distance 
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graphs. A dense distance graph for a cluster C is defined to be a complete graph over the border 
vertices of C where edge lengths correspond to shortest path distances in C. In order to compute 
dense distance graphs for all clusters we use Klein's algorithm [11], who have shown that after 
O(nlogn) preprocessing time any distance from the external face can be computed in O(logra) 
time. The proof of the following lemma is given in Appendix [Bl 

Lemma 8. Given an r -partition V we can compute a dense distance graph for all clusters in V in 
0{n log r) time. 

The dense distance graphs can be used to speed up shortest path computations using Dijkstra's 
algorithm. It was shown by Fakcharoenphol and Rao ([2], Section 3.2.2) that a Dijkstra-like algo- 
rithm can be executed on a dense distance graph for a cluster P in O(\0P\ log 2 \P\) time. Having 
constructed the dense distance graphs, we can run Dijkstra in time almost proportional to the 
number of vertices (rather than to the number of edges, as in standard Dijkstra). We use this 
algorithm in graphs composed of dense distance graphs and a subset E' of edges of the original 
graph G = (V, E): 

Corollary 9. Dijkstra can implemented in 0{\E'\ log \V\ + ^^dG^lo^ \dGi\) time on a graph 
composed of a set of dense distance graphs Gi and a set of edges E' over the vertex set V . 

Proof. In order to achieve this running time we use Fakcharoenphol and Rao [2] data-structure for 
each Gi. Moreover, minimum distance vertices from each Gi and all endpoints of edges in E' are 
kept in a global heap. □ 

In general, clusters may contain holes although in typical cases, e.g., in grid graphs, the obtained 
clusters are holeless. In this section, in order to introduce the main ideas behind our algorithm, 
we restrict ourselves to holeless r-partitions, i.e., r-partitions where each cluster contains one hole 
(external face). We will show in Appendix ID1 how to handle the general case. 

Assume that we have computed the dense distance graphs for all clusters in a given r-partition 
V of the dual graph Gr>. Recall that in our min-cut algorithm we are free to choose any path n 
from f s to ft in Go, as long as tt is non-crossing. We choose n to minimize the number of clusters 
it crosses and to maximally use the dense representation of clusters. 

We define a skeleton graph Gp = (dV, Ep) to be a graph over the set of border vertices in V. 
The edge set E-p is composed of infinite length edges connecting consecutive (in the order on the 
hole) border vertices on each hole. By our holeless assumption, all border vertices in each cluster 
lie on the external face of the cluster, so the graph G-p is connected. We define a patched graph 
to be G = Go U G-p. Note that this graph is still planar and the shortest distances do not change 
after adding infinite length edges. 

Define G to be the graph composed of: (1) two clusters P s and Pt that include f s and ft 

respectively; (2) the dense distance graphs (represented by square matrices) for all other clusters; 

s t 

(3) the skeleton graph Gp (see Figure [6|). Note that G ' contains: 0(2r + jy/r) = 0(r + -^=) 

vertices; 0(j) dense distance graphs each over 0(^/r ) vertices; at most 3r edges of Gd from P s 
and P t ; additional O(^) edges from the skeleton graph Gp. Using Corollary [9] to run Dijkstra's 
algorithm we get the following. 

Corollary 10. The shortest paths in G S,t can be computed in 0((r + -^)log 2 n) time. 
3.1 Recursive Division 

Let b s and bt be any border vertices in clusters P s and Pt respectively. We define it to be com- 
posed of: a simple path from f s to b s in P s ; a simple path from b s to bt in Gp; and a simple path 
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from bt to ft in P%. Note that the construction of the graph G ' and consequently of takes 
0(n log r + ^= log n) time by Lemmas [7] and [H 

Let Q be some /j-cut-cycle in Cp. After finding Cj, we need to recurse on the graphs Gofl 
int(Ci) and Gofl ext{Ci). These graphs cannot be computed explicitly. However, we will show how 
to determine G := Gp> n mi(Cj) ' and G := C?d n ext{Ci) ' using G ' . Let P be a cluster in 
the partition other then P s or P^, and let Gp be its dense distance graph. For a set of vertices X we 
define GpCiX to be the dense distance graph of P(~)X . The /j-cut-cycle was found using G ' so parts 
of Ci that pass through P correspond to shortest paths. Hence, the shortest paths in Pr\int(Ci) and 
PC\ext{Ci) cannot cross the cycle Cj. As a result, distances in GpC\int{Ci) and Gpr\ext(Ci) between 
border vertices of P that are not separated by Cj are the same as in Gp. On the other hand, for 
border vertices that are separated by Cj the distances are infinite (see Figure [7]). We define G ' OX 
to be the graph obtained by taking GpDX for every cluster in G. Note that we have the following. 

Corollary 11. G <1 = (?'* n int(Ci) and = Zf '* n ext{d). 

Using these equalities the construction of G <1 and G l> takes time proportional to the size of G s,t 
only. Given a path ir = fx, . . . , fk, the recursive algorithm for computing min-cuts works as follows. 

S t 

First, remove from G ' vertices with degree 2 by merging the two incident edges, and find an /j-cut- 
cycle C{ for i = [^\ ■ Next, construct graphs G and G , and recursively find a minimum cut C<j 

in G <1 and a minimum cut Cj> in (T > . Finally, return the smallest of the three cuts C<j, Ci and Cj>. 

To achieve our promised bounds, we need to show that the above algorithm works in sublinear 
time. In order to show that, we only need to prove that the total size of the graphs on each level 

S t 

in the recursion tree is small. Take a graph G ' and for each cluster P that contains more then 
one border vertex add a new vertex vp and replace each dense distance graph by a star graph with 

S t 

in a center. We call the resulting graph the contracted graph and denote it by G c ' . Obviously, 
the contracted graph has more vertices than the original graph: thus, we can bound the number 
of vertices in the original graph by considering only contracted graphs. The proof of this lemma is 
included in Appendix O 

Lemma 12. The total number of vertices in contracted graphs on each level in the recursion tree 
is 0(r + f r ). 

By the above lemma and Corollary [91 running Dijkstra for each level takes 0((r + -^=)log 2 n) 
time in total. On each level the length of the path ir is halved, so there are at most logra recursion 
levels. Hence we obtain the following theorem. 

Theorem 13. Let G be a flow network with the source s and the sink t. If Go allows holeless 
r -partition, then the minimum cut between s and t can be computed in 0(n log r + (r + log 3 n) 

time. By setting r = log 8 n we obtain an O(nloglogn) time algorithm. 

Theorem 1131 holds also for general r-partitions within the same 0(n log log n) time bound. The 
modifications needed to make the above algorithm work in the general case are presented in Ap- 
pendix |Dj 

4 Computing Maximum Flows 

The standard ways of computing maximum flow in near linear time assume that we already have 
computed its flow value /. It is given by the min-cut capacity and as we already know it can be 
computed in 0(n log log n) time. We will use the approach proposed by Hassin and Johnson [7], 
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but adopted to our case as we use a different family of cuts Ci. As argued in the following, their 
approach uses only very basic properties of the cut-cycles and can be directly applied to our case. 
Moreover, for the sake of brevity, we will assume that that we are given a holeless r-partition. The 
general case can be handled using ideas presented in Appendix [Dl and it will be included in the 
full version of this paper. 

Let us define the graph to be the graph obtained from G' n by adding directed edges of length 
— / from fi to f- for all 1 < i < k. After Lemma 4.1 in Miller and Naor [12] we know the following. 

Corollary 14. The graph does not contain negative length cycles. 

Let r be such that C r is the shortest of all Ci for 1 < i < k. The above corollary assures that 
distance 5(v) from f' r to a vertex v is well defined in G^. The next lemma follows by Theorem 1 
in [7J or Section 5.1 in [12] , 

Lemma 15. Let e = (u,v) be the edge in G and let er> = (/«, f v ) be the corresponding edge in Gd- 
The face f u is defined to lie to the left when going from u to v. The function f(u, v) := 5(f u ) — 5(f v ) 
defines the maximum flow in G. 

By this lemma, in order to construct the flow function we only need to compute distances from 
f' r in ,r. The cycles Gj for i = 1, . . . , k divide into k + 1 subnetworks C^Oi • • • ; ~G* k where for 
1 < i < k — 1, is the subnetwork bounded by and including Ci and Cj+x- The next lemma is a 
restatement of Lemma 2 from [7J . 

Lemma 16. Let v be a vertex in iVj. Then if i < r then there exists a shortest (f^.,v)-path which is 
contained in U}=o^j- Similarly, if i > r then there exists a shortest {f' r ,v)-path which is contained 
in [f j=r Gj. 

Lemma [TBI implies that the computation of 5(v) for v G U^=o can be restricted to |J^=o 

„ nl , We can r e strict th e coition for „ 6 U*„ % ^ * ^ Won. 

Similarly to [7J let us define a normal path to be a simple {f' T , v )-path p(v) = p r ■ . . . ■ p q ■ . . . ■ P2q-i 
such that, for j = r, . . . , q, subpath pj is in Gj, and, for j = q + 1, . . . , 2q — i, subpath pj is in 

G^2(7-j and uses no edges of negative length. We require as well that q is minimal. Hassin and 
Johnson [7J have shown that distances from f' r in an n- vertex graph can be computed in O(nlogn) 
time when for each vertex v there exists a shortest {f' r , f)-path that is normal (Theorem 2 in [7J). 
Their computation is based on Dijkstra's algorithm and can be sped up using dense distance graphs 
and Corollary [9j 

Lemma 17. Distances from f r in the graph ~dt n can be computed in 0(nlogr + + r) log 2 n) 
time when for each vertex v there exists a shortest {f' r ,v)-path that is normal. 

Proof. The algorithm works in two phases. First, we take and substitute each cluster P with 
its dense distance graph. The resulting graph has 0(-^) vertices. The Dijkstra-like computation 

from [7J can be executed in 0((-^ + r) log 2 n) time on this graph. In this way we obtain distances 
5(v) for border nodes v in all clusters in Gd. Second, for each cluster separately we run Hassin and 
Johnson's computation starting from border vertices only. The second phase works in 0(r log r) 
time for each cluster, which yields 0(n log r) time in total. □ 

In order to use the above lemma we only need the following result, which can be proven in the 
same way as Lemma 3 from [7J . The proof in [7j uses only the fact that subpaths of cut cycles Ci 
are shortest paths and this holds in our case as well. 
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Lemma 18. For any i = r, . . . , k, for every vertex v, in £r , there exists a normal shortest path p(v). 

Proof. Assume that a shortest path p(v) intersect some cycle Cj. If the previous intersected cycle 
was Ci as well, then we can either short cut p(v) using the part of Ci or short cut Ci using part 
of p(v) (see Figure [8]). This contradicts either the minimality of p(v) or the minimality of C{. 
Moreover, after p{v) leaves Ci to go into ~G*i-\ it cannot use a negative edge any more, as otherwise 
it would cross itself. If it crosses itself then the resulting cycle could be removed as it cannot have 
a negative weight. □ 

Combining Lemma [T5"l Lemma [T71 and Lemma ITHI we obtain the main result of this section. 

Theorem 19. The maximum flow in an undirected planar graph can be computed in 0(n log r + 
(■^ + r) log 2 n) time. By setting r = log 8 n we obtain an O(nloglogn) time algorithm. 

5 Dynamic Shortest Paths and Max Flows in Planar Graphs 

Most of the ideas presented in this section are not entirely new, but nevertheless combined together 
they are able to simplify and improve previously known approaches. Our dynamic algorithm builds 
upon the r-partition introduced in Section [3l We first show how to maintain a planar graph Gd 
with positive edge weights under an intermixed sequence of the following operations 0: 

insert£)(x, y, c) add to Gd an edge of length c between vertex x and vertex y, provided that 

the new edge preserves the planarity of Gd; 
deleteci(x,y) delete from Gd the edge between vertex x and vertex y; 

shortest_path[)(x,y) return a shortest path in Gd from vertex x to vertex y. 

In our dynamic algorithm we maintain the r-partition of Gd together with dense distance graphs 
for all clusters in the partition. This information will be recomputed every \/r operations. We now 
show how to handle the different operations. We start with operation insert: let (x,y) be the edge 
to be inserted, and let P x (respectively P y ) be the cluster containing vertex x (respectively vertex 
y). If x and y are not already border vertices in clusters P x and P y , we make them border vertices 
in both clusters and add edge (x,y) arbitrarily either to cluster P x or to cluster P y . Next, we 
recompute the dense distance graphs of clusters P x and P y , as explained in Lemma [HI This requires 
overall time 0(r log r). Note that an insert operation may increase by a constant the number of 
border vertices in at most two clusters and adds an edge to one cluster: since the partition into 
clusters is recomputed every y/r operations, this will guarantee that throughout the sequence of 
operations each cluster will always have at most 0(r) edges and 0(y/r) border vertices. To delete 
edge (x,y), we remove this edge from the only cluster containing it, and recompute the dense dis- 
tance graph of this cluster. Once again, this can be carried out in time 0(r log r). Amortizing the 

nlogr+-?= logn 

initialization over \Jr operations yields 0{ -^=f- V r logr) time per update. 

— s t 

In order to answer a shortest_path(x , y) query, we construct the graph G ' (as defined in Sec- 

s i 

tion[3j). Note that the distance from s to t in Gd and G ' are equal. Hence, by Corollary 1101 the 
shortest path from s to t can be computed in 0((r + ^j=) log 2 n) time. 

In order to minimize the update time we set r = re 2 / 3 and obtain the following theorem. 



1 We have chosen to work with the dual graph Gd for consistency with the remaining parts of the paper. 
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Lemma 20. Given a planar graph G with positive edge weights, we can insert edges (allowing 
changes of the embedding), delete edges and report shortest paths between any pair of vertices in 
0(n 2 / 3 i Q g2 n ^ amor ti ze d time per operation. H 

We recall that we can check whether a new edge insertion violates planarity within the bounds 
of Lemma [20l indeed the algorithm of Eppstein et al. [T] is able to maintain a planar graph subject 
to edge insertions and deletions that preserve planarity, and allow to test whether a new edge 
would violate planarity in time 0(n 1//2 ) per operation. Finally, we observe that not only we have 
improved slightly the running time over the 0(n 2 / 3 log 7 / 3 n) time algorithm by Fakcharoenphol and 
Rao [2], but our algorithm is also more general. Indeed, the algorithm in [2] can only handle edge 
cost updates and it is not allowed to change the embedding of the graph. On the contrary, our 
algorithm can handle the full repertoire of updates (i.e., edge insertions and deletions) and it allows 
the embedding of the graph to change throughout the sequence of updates. 

We now turn our attention to dynamic max-flow problems. In particular, given a planar graph 
G = (V, E) we wish to perform an intermixed sequence of the following operations: 

insert(x, y, c) add to G an edge of capacity c between vertex x and vertex y, provided that 

the embedding of G does not change; 
delete{x, y) delete from G the edge between vertex x and vertex y; 

max-flow(s,t) return the value of the maximum flow from vertex s to vertex t in G. 

Note that insert operations are now more restricted than before, as they are not allowed to 
change the embedding of the graph. To answer max-flow queries in the primal graph G we need to 
maintain dynamically information about the distances in the dual graph Gd, with the bounds re- 
ported in Lemma [201 Unfortunately, things are more involved as a single edge change in the primal 
graph G may cause more complicated changes in the dual graph Gd- In particular, inserting a new 
edge into the primal graph G results in splitting into two a vertex in the dual graph Gd, whereas 
deleting an edge in the primal graph G implies joining two vertices of Gd into one. However, as 
edges are inserted into or deleted from the primal graph, vertices in the dual graph are split or 
joined according to the embedding of their edges. To handle efficiently vertex splits and joins in 
the dual graph, we do the following. Let / be a vertex of degree d in the dual graph: we maintain 
vertex / as a cycle Cfoid edges, each of cost 0. The actual edges originally incident to /, are made 
incident to one of the vertices in Ct in the exact order given by the embedding. It is now easy to 
see that in order to join two vertices fi and /2, we need to cut their cycles Cf 1 and Cf 2 , and join 
them accordingly. This can be implemented in a constant number of edge insertions and deletions. 
Similarly, we can support vertex splitting with a constant number of edge insertions and deletions. 
Additionally, for each cluster P, for each pair h, hi of holes we need to compute the dense distance 
graphs of P-w hh , and the minimum-cuts between bh and by. However, following Corollary 1241 this 
does not increase the running time of our dynamic algorithm. Note that this information is enough 
to construct the graphs G s ' and G% in 0(r + ^) time. Moreover, by Theorem 1131 and Theorem [26l 

the min-cut algorithm can be executed on these graphs in 0((r + log 3 n) time. Setting r = n 2 / 3 
this yields immediately the main result of this section. 

Lemma 21. Given a planar graph G with positive capacities, each operation insert, delete and 
max_flow can be implemented in 0(n 2//3 log 3 n) amortized time. H 



2 The same worst-case time bounds can be obtained using a global rebuilding technique. 

3 Again, the same worst-case time bounds can be obtained using a global rebuilding technique. 
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Figure 1: An embedded planar graph G and its dual graph Go- Vertices of G are shown as circles 
and edges of G are solid. Vertices of Go are shown as gray squares, and edges of Go are dashed, t 
and t are two vertices in G, and f s and ft are arbitrary inner faces incident respectively to s and t 
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Figure 2: (a) The dual graph Gd of Figure HJ the path ir from / s to ft is shown in bold, (b) 
The dual graph Gd embedded so that the path ir is a horizontal line, (c) The graph G' n . (d) The 
graph G'l = G' w , embedded so that the face defined by ir and ir' is the outer face, (e) The graph G w 
obtained after identifying the path n' (respectively n) in G" with the path it (respectively ir 1 ) in G' n . 
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(e) - - 

(f) 

Figure 3: On the proof of Lemma [3j (a) The /j-cut-cycle C crosses the path tt in in the graph 
Gd- (b) The corresponding /^-separating path p in G w . (c) The /j-cut-cycle C touches from 
above the path 7r in /, in the graph Gd- (d) The corresponding /j-separating path p in G n . (e) 
The /j-cut-cycle C touches from below the path tt in /j in the graph Gi> (f) The corresponding 
/i-separating path p in G^. 




Figure 4: On the proof of Lemma [H 



14 



(a) (b) 

Figure 5: Dealing with non-simple paths, (a) A non-crossing path ir of Gp> that touches itself at 
face v = V{. (b) The vertex splitting transformation on face v. 




— s t 

Figure 6: The graph G ' containing clusters P s and Pt, the dense distance graphs for all other 
clusters, and the skeleton graph G-p. The path ir from f s to ft is shown in bold. 




C, 



Figure 7: Computing Gp D int(Ci). It two border vertices u and v are not separated by the 
/j-cut-cycle Cj, then their distance in Gp n int(Ci) is equal to their distance in Gp. On the 
other hand, if border vertices v and w are separated by the /j-cut-cycle Cj, then their distance in 
Gp n int(Ci) is infinite. 
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A Computation of r-Partition: Proof of Lemma [7] 



In order to prove Lemma [7] we will combine Frederickson's [5] technique with the following impli- 
cation of Section 3.1 and Section 5.1 from [2]. 

Corollary 22. An r-partition of an n-node planar graph that already contains 0(^j=) border ver- 
tices can be computed in 0(n log n) time. 

Proof. Actually, Fakcharoenphol and Rao [2] construct in the above time bound a recursive decom- 
position and in order to get an r-partition we only need to run their algorithm until the size of the 
clusters drops below 0(r). 

Moreover, the recursive decomposition of Fakcharoenphol and Rao [2] assumes that border ver- 
tices are present in the decomposed graphs or clusters. Moreover, each time a cluster is split into 
smaller clusters, the border vertices are split into asymptotically equal parts. The 0(-^=) border 
vertices we start with will be distributed equally into O(^) clusters. Hence, each cluster in the 
obtained partition will have additional 0(-p= x ~) = 0{y/r) border vertices. □ 

Now in order to the find the r-partition quickly we use the following algorithm: 

• Generate a spanning tree T of the graph; 

• Find connected subtrees of T containing 0{y/r) vertices using a procedure from [S]; 

• Contract the graph on these subsets, to obtain a simple planar graph G s with 0(-^=) vertices; 

• Using Corollary 1221 find an r-division in G s with 0(^372) clusters of size O(r); 

• Expand G s back to G. In G there are 0{n/r) clusters V\ of size 0{\/r) resulting from boundary 
vertices in G s , and 0(^372) clusters Vi of size 0(r 3 / 2 ) resulting from the interior vertices of G s \ 

• Apply Corollary [22] to find an r-division for clusters in Vi taking into account the border vertices 
already present in the clusters. 

It is easy to see that the above procedure requires 0(n log r + -^j= logn) time. We only need to show 
that the result is a valid r-division. The clusters V\ correspond exactly to the connected subtrees 
of T of size 0(y/r). They cannot have more then 0(y/r) border vertices and all the border vertices 
lie on the external face, i.e., they contain one hole. Hence, clusters in V\ satisfy the properties of 
an r-partition. 

Consider the clusters in Vi before we have applied Corollary [22j Each cluster is obtained by 
expanding a cluster in G s that had \fr border vertices. Consider the process of expanding to a 
subtree T a border vertex bh lying on a hole h in a cluster P. Let h' be the hole obtained from h 
by removing b^ from P. The edges of T were not present in the cluster so the process of expanding 
T can be seen as gluing T at the side of U . Note that by the connectivity of T no new hole is 
created and only the hole h' becomes "smaller". Moreover, not all of the vertices of T become 
border vertices of the expanded piece, as some of them do not lie on the side of the new hole (see 
Figure [9j). Nevertheless as T contains 0(y / r) vertices at most 0(y / r ) vertices can become border 
vertices. Hence, in total the pieces in Vi will have O(r) border vertices. This number satisfies the 
assumption of Corollary [221 so pieces obtained by using it satisfy the assumption of r-partition. 
This completes the proof of Lemma [7J 

B Computation of Dense Distance Graphs: Proof of Lemma [8] 

Proof. In order to compute dense distance graphs we use the following result by Klein |llj . 

Theorem 23 (Klein |llj). Given an n-node planar graph with non-negative edge lengths, it takes 
0(n\ogn) time to construct a data structure that supports queries of the following form in O(logn) 
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Figure 9: a) Before expanding the border vertex bh in hole h; b) after expanding bh to T no new 
hole is created. 



time: given a destination vertex t on the boundary of the external face, and given a start vertex s 
anywhere, find the distance from s to t. 

A cluster C in the r-partition has r vertices, 0(y/r) border vertices and a constant number of 
holes. Border vertices in C lie on one of the holes. In order to compute distances from all border 
vertices to a given hole H we will apply Theorem 1231 We simply find an embedding of the graph 
such that H becomes the external face, and query the distances from all border vertices to H. 
There are a constant number of holes in each cluster, and so Klein's data structure will be used 
a constant number of times for each cluster. On the other hand, there are 0(r) pairs of border 
vertices in each cluster, so we will ask 0(r) queries. The time needed to process each cluster is 
hence 0(r log r), which gives ^ x 0(r log r) = 0(n log r) time in total. □ 



C Size of Contracted Graphs: Proof of Lemma U2 



Proof. At the top level of the recursion we have one graph G*' that has 0(r + -^j=) vertices and 
0{r + ^j=) faces. Now consider the case when the graph G s ' 1 is split along Cj into G <% and G l> . 

— s t 

Observe that the cycle Cj can be traced in G c ' when we replace each edge in dense distance graphs 
for a cluster P by a two-edge path going through the center vertex vp. In such a case we have 

^2 S t Z~> S t 

G c = G c ' n int(Ci) and G c = G c ' n ext(Ci). In other words each time we recurse we split the 

— s t 

graph G c ' along some cycle and then remove degree-two vertices. Note that after splitting the graph 

— s t 

G c ' along the cycle Cj, the number of faces in the union of the resulting graphs increases by exactly 
one (see Figure [T0|) . There are at most 0(r + -^) recursive calls, so the union all contracted graphs 
on a given level in the recursive tree has 0(r+-^) faces as well. Moreover, the vertices in this union 
graph have degree at least three, because all degree-two vertices are removed. Let v, e and / be the 
number of vertices, edges and respectively faces in the union graph. Now, by Euler's formula the 
claim of the lemma follows: v = 3v-2v = 3v-2(2+e-f) < 2e-2(2+e-/) < 2/-4 = 0(r+-^). □ 



D General Clusters 

When holes are present in the r-partition the skeleton graph G-p is not connected. Hence, the path 
connecting b s and bt cannot use border vertices only. We need to modify the algorithm to allow 
such paths. We do this as follows. 
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(a) 




(b) (c) 
s i 

Figure 10: (a) The graph G c ' and the cycle Cj. Cluster boundaries are shown in bold, and 

— i> — s t 

the cycle Cj is shown with dashed edges, (b) The graph G c = G c ' n ext{Ci). (c) The graph 
G^ 1 = G^.' n int(Ci). Note that the union of G^ and G^T* contains one more face than G s c ' . 
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For each hole h in P, we fix a border vertex bh- For each pair of holes h, h! in P, we fix a path 
TTh,h' that starts from bh and ends in bh>, goes through bh" for all holes h" in P, and for all bh" on 
the path walks around the hole h" passing through all its border vertices (see Figure [TT]) . These 
paths are used to do some additional preprocessing for each cluster P in the partition. For each 
pair of holes h, h! in P we compute dense distance graph for P 7Th , , and find the minimum cut Ch,h' 
between bh and bh' in P. 



h 









h' } 


( h" ■ 



Figure 11: The path ~Kh,h'- 



Corollary 24. The additional preprocessing takes 0(n log r) time. 

Proof. For each cluster P and each pair of holes the dense distance graph can be computed in the 
same manner as in Lemma On the other hand, minimum-cuts can be found in 0{r log r) time 
using the algorithm by Henzinger et al. [8]. Hence, over all clusters we need a total of 0{n log r) 
time. □ 

Now in order to connect f s and ft we will use paths TThh' whenever we need to pass between 
two different holes h, h! in a piece P. Let V n be the set of all such pieces on tt. The resulting path 
is no longer simple, but it will be non-crossing. As shown in Section \2.2\ our min-cut algorithm can 
be executed on non-crossing paths as well. We can make the following observation (see Figure fT2j) . 

Corollary 25. A minimum s-t cut C either contains a vertex in dir or is fully contained in one 
of the pieces . 

Proof. If the cycle C contains a vertex in chr, we are done. Assume that it does not contain any 
vertex from dir. In order to be a cut, it has to cross path n and it can do it so in one of nh,h'- 
Then it has to be fully contained in the corresponding P as by the construction of h> or border 
vertices of P lie on ir. □ 

Using Corollary [25l we can find the minimum cut in two phases. First, let Cj be the smallest 
of the cuts Ch : h' in 'P-k for bh,b' h E n. Second, run the algorithm from previous section on a path 

— s t 

dir in G to find a cut Cf,. Finally, return the smallest of the cuts Cj and CV The running time 
of the above algorithm is the same as in Theorem [13j 

Theorem 26. Let G be a flow network with source s and sink t. The minimum cut between s and 
t in G can be computed in 0(n log r + (r + -^=) log 3 n) time. By setting r = log 8 n, we obtain an 
O(nloglogn) time algorithm. 
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Figure 12: On illustrating Corollary [25j A cut can be either fully contained in one of the pieces 
V n (such as cut Cj) or it must contain one of the border vertices in dir (such as cut Cj). 
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